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Abstract. We report the release of the Gauss-Legendre Sky Pixelization (GLESP) software package version 
1.0. In this report we present the main features and functions for processing and manipulation of sky signals. 
Features for CMB polarization is underway and to be incorporated in a future release. Interested readers can visit 
http://www.glesp.nbi.dk/ and register for receiving the package. 
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1. Introduction 

The CMB temperature fluctuations on a sphere can be 
written as a sum over spherical harmonics 



(1) 
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where \ai m \ are the moduli of the spherical harmonic co- 
efficients of the expansion and 4>i m are the phases defined 
in [0, 2tt]. The angular power spectrum is defined as 
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where T is the CMB temperature 2.726 K. 

In this article, we report the first release of 
the Gauss-Legendre Sky Pixelization (GLESP) soft- 
ware package for the CMB analysis. The GLESP 
scheme UDoroshkevich et al. 2005(1 is based on the Gauss- 
Legendre polynomial zeros, allowing one to perform strict 
orthogonal spherical harmonic expansion of the CMB 
temperature anisotropics. The GLESP package has been 
tested on the following operating systems: RedHat Linux, 
SGI Irix and Solaris. Features for CMB polarization are 
under preparation. 



2. Features in the GLESP package vl.O 

2.1. Basic features for CMB analysis 

The basic features of the GLESP package includes spher- 
ical harmonic analysis of a map, map synthesis from a 
given set of a\ m coefficients, and simulation of a Gaussian 
map from a given angular power spectrum Ci , or Di . 

cl2map 

cl2map synthesizes CMB temperature maps from 
spherical harmonic coefficients a; m . Inversely, it can 
decompose an input CMB temperature map into its 
spherical harmonic coefficients ai m and/or power spec- 
trum. If the input is a power spectrum C/, or Di, 
cl2map generates a Gaussian map of CMB anisotropics 
by assigning the ai m mutually independent Gaussian- 
distributed real and imaginary parts, both with a stan- 
dard deviation y/Ci/2 at each I harmonic, 
rsalm 

rsalm convolves or deconvolves a map with a symmet- 
ric beam via 



— Bi aw, 



(4) 



where B\ is the expansion in Legendre polynomials of 
a symmetric beam. The symmetric beam can be in- 
put by reading in an ASCII file containing the pro- 
file of B\ or by providing the F WHM for a symmetric 
Gaussian beam, where Bi = exp[— a 2 l(l + l)/2] and 
a = FWHM/V8hi2. 
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For deconvolution, a simple Tikhonov regularization 
scheme is applied via 



air, 



Bi+a 



(5) 



where a is the regularization parameter. 
f2flg 

f 2f ig converts a GLESP map in FITS binary format 
into GIF (Graphic Interchange Format). f2f ig offers 
users to put markers on the output figure. The mark- 
ers can be input with one of the following 3 coordinate 
systems: in 9 and <f> of spherical polar coordinate sys- 
tem, in Galactic latitude b and longitude I of Galactic 
coordinate system and in right ascension (R.A.) and 
declination (Dec.) of equatorial coordinate system. 
alm2dl 

alm2dl reads a/ m and calculates the angular power 
spectrum Ci or D[. It also provides the options in 
listing the amplitudes \ai m \ or phases <pi m defined in 
[0,2n]. 

2.2. Advanced features: map and ai rn manipulation 

The package offers several advanced features that allow 
the users to manipulate pixelized maps or sets of et/ m . 

difmap 

dif map performs linear combinations of maps, map ro- 
tation, conversion of a map from Galactic to equatorial 
coordinate base. It can display simple map informa- 
tion such as pixel numbers and resolution, histogram, 
mean and variance, maximum and minimum temper- 
ature values and positions in (9, <j>) . 
difalm 

dif aim performs arithmetic operations over et; m such 
as addition, subtraction, multiplication and division 
by a factor. It also allows to swap the amplitudes and 
phases between two sets of a/ m . It can multiply /divide 
o-im by given vectors v(l). 
mapcut 

mapcut modifies a map in various ways. It can sin- 
gle out an area (either circular or rectangular) and 
set temperature values to zero in the area or outside 
the area. Furthermore, it can systematically remove 
or withhold several circular (or rectangular) areas by 
reading in a file containing those areas' positions and 
radii (or two corners). It can dissect a rectangular area 
from a map and output to a file in FITS format, 
mapcut can also display temperature values of a 
GLESP map at given positions in one of the 3 coordi- 
nate systems: spherical polar, Galactic, and equatorial, 
mappat 

mappat performs the conversion between ASCII and 
FITS binary formats for a GLESP map. The derivative 
of this is that it allows one to input any patterns such 
as point source positions in ASCII and produce their 
whole-sky distribution. 



2.3. Conversion between HEALPix and GLESP 

The package also provides map conversion between 
HEALPix flGorski, Hivon fc Wandelt 1999D and GLESP. 

f2map 

f2map performs repixelization of a GLESP map to 
HEALPix or to equidistant grid FITS Basic format. 
In the process of conversion to HEALPix, f 2map can 
perform repixelization cither using spline interpolation 
algorithm or direct mapping, 
cmap 

cmap performs repixelization of a HEALPix map to 
GLESP using spline interpolation for re-pixelization. 
On conversion to GLESP, the resolution is determined 
by the N S id e value in HEALPix, otherwise, the resolu- 
tion can be set according to user's input, 
cmap can also display temperature values of a 
HEALPix map, or output the HEALPix map in ASCII 
format in 3 columns: 9, (j) and temperature, cmap 
also reads the 2nd field of Binary Table Extension in 
HEALPix maps, 
ntot 

ntot calculates the total pixel number and size for 
both GLESP and HEALPix schemes. By giving a res- 
olution, ntot calculates the parameters of the nearest 
allowed resolution for both schemes. 

3. Objectives 

3.1. Operation on a single map 

Decomposing a map for its a; TO /C/ cl2map 
Simulating a Gaussian map from given C; cl2map 
Displaying the map T max /T m i ni map resolution 

mean/variance difmap 
Map rotation difmap 
Arithmetic operations difmap 
Producing the histogram difmap 
Converting the map coordinate system to 

equatorial coordinates difmap 
Casting a mask on a map difmap 
Plotting the map and putting markers on the 

figure f2fig 
Displaying the temperature value at a given 

coordinate mapcut 
Truncating the temperature fluctuation to 

a given temperature range mapcut 
Removing or keeping specified areas from 

the map mapcut 
Map convolution and deconvolution rsalm 



3.2. operation on a single set of a\ m 

Calculating the amplitude |aj m | alm2dl 

Calculating the phases <pi m alm2dl 

Synthesizing a\ m for the map cl2map 

Arithmetic operations on a; m difalm 
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Shifting phases dif aim 

3.3. Operations over several maps 

Linear combinations over maps dif map 

Cross correlation of 2 maps dif map 

Multiplication of 2 maps dif map 

3.4. Operation over several sets of a\ m 

Linear combinations of ai rn dif aim 

Swapping phases and amplitudes of 2 ai m dif aim 

3.5. Conversions 

Converting GLESP FITS to HEALPix FITS cmap 

Converting HEALPix FITS to GLESP FITS f 2map 

Converting FITS binary to ASCII mappat 

Converting/creating FITS binary from ASCII mappat 

Information on pixel number and size ntot 



4. Acknowledgment of using GLESP software 
package 

Those who use the GLESP package for producing re- 
sults should acknowledge in the publication. The refer- 
ence should include the paper introducing the GLESP con- 
cept : Doroshkevich et al. International Journal of Modern 
Physics D, Vol 14, No.*, *** (2005). 

5. Problems and Suggestions 

For any technical issues, bugs, problems or suggestions, 
please send email to glesp@nbi.dk. 
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